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^ ■ Abstract 

Exclusive processes at large momentum transfer Q factor into perturba- 
tively calculable short-distance parts and long-distance hadronic wave func- 
tions. Usually, only contributions from the leading Fock states have to be 
included to leading order in 1/Q. We show that for exclusive decays of 
^ P-wave quarkonia the contribution from the next-higher Fock state |QQg) 

contributes at the same order in 1/Q. We investigate how the constituent 
gluon attaches to the hard process in order to form colour-singlet final-state 
hadrons and argue that a single additional long-distance factor is sufficient 
to parametrize the size of its contribution. Incorporating transverse degrees 
of freedom and Sudakov factors, our results are perturbatively stable in the 
sense that soft phase-space contributions are largely suppressed. Explicit 
calculations yield good agreement with data on Xcj decays into pairs of 
pions, kaons, and etas. We also comment on J/i[) decays into two pions. 
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1 Introduction 



Exclusive reactions at large momentum transfer Q can be calculated in perturbative QCD 
(pQCD) owing to a factorization theorem |TJ], which separates the short-distance physics 
of the partonic subreactions at the scale Q from the longer-distance physics associated 
with the binding of the partons inside the hadrons. The full amplitude is given as a sum of 
terms, where each term factors into two parts, a hard-scattering amplitude T# , calculable 
in perturbative QCD, and wave functions ip(xi, kj_i) for each hadron H. The importance of 
the various terms depends on their scaling with 1/Q. For the contribution with the weakest 
fall-off with Q (leading- twist contribution), the amplitude T# describes the scattering of 
clusters of collinear partons from the hadron and is given by valence-parton scatterings 
only. Hence the only non-perturbative input required are the distribution amplitudes 
4>H{%i, Q) for finding valence quarks in the hadron, each carrying some fraction Xj of the 
hadron's momentum. The distribution amplitudes represent wave functions integrated 
over transverse momentum up to a factorization scale fip of order Q. 

Corrections to this standard hard-scattering approach (sHSA) (lowest-order in a s (Q) 
calculations in the collinear approximation using valence Fock states only) can be of 
perturbative origin (oc a s (Q)) or power-like (oc 1/Q). The latter may be classified as 
follows: 

1. Corrections arising from the overlap of the soft wave functions. 

2. Corrections associated with the transverse momentum of the partons inside the 
hadrons. 

3. Corrections from higher Fock states. 

In order to show that predictions of the sHSA are reliable, two conditions have to be 
met. First, one has to prove that both the 0(a s (Q)) corrections and the 1/Q corrections 
are small. Second, one has to ensure that the perturbative calculation of Th does not 
receive major contributions from phase space regions, where the virtualities of the internal 
partons become soft and perturbation theory is not applicable @. 

In recent years there has been progress in our understanding which of the above- 
mentioned corrections are important. On the one hand, the HSA has been modified 
through the incorporation of tranverse degrees of freedom and gluonic radiative corrections 
in form of a Sudakov factor ||. Within this modified hard-scattering approach (mHSA), 
the perturbative contribution to exclusive observables can be calculated self-consistently, 
because the Sudakov factor suppresses contributions from soft phase-space regions. In 
particular, the mHSA provides a description of the pion-photon transition form factor, 
which is both reliable and in agreement with data, already at momentum transfers as low 
as lGeV [§§,§]. 

On the other hand, the overlap of the initial- and final-state pion wave functions, 
representing a soft contribution of higher-twist type, is still sizeable in the few GeV 
region in the calculation of the pion elastic form factor and just suffices to fill the gap 
between the perturbative prediction and the experimental data (7], |J. Since the overlap 
contribution, while commonly neglected, is an essential ingredient of the HSA [|I|], the 
smallness of the perturbative contribution to the pion form factor should not be regarded 
as a failure but rather as consistent within the entire approach. 

In a recent paper || we have shown that the hierarchy of the 1/Q expansion for 
exclusive decays of heavy quarkonia is different from the canonical one: for P-wave states, 
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the contribution from the next-higher Fock state is not suppressed by additional powers 
of 1/Q compared to the contribution from the leading Fock state. (Q is of the order 
of the heavy-quark mass m c .) This is analogous to the case of inclusive decays of heavy 
quarkonia, although there the long-distance matrix elements are organized into a hierarchy 
according to their scaling with v, the typical velocity of the heavy quark in the quarkonium 

I- 

For the P-wave charmonium state Xcj, the Fock-state expansion starts as 

\XcA J ++ )) = 0(1) MPj)) + 0(v) |cc 8 ( 3 ^) g) + 0(v 2 ) , (1.1) 

where the subscript c specifies whether the cc c is in a colour-singlet (c = 1) or colour- 
octet (c = 8) state. The angular-momentum state of the cc pair is denoted by the 
spectroscopic notation 2S+1 Lj. Contrary to the case of S"-wave decays, two 4-fermion 
operators contribute to the decay rate of P-wave states into light hadrons at leading 
order in v 

T[ XcJ - LH] = \ (Xcj|d( 3 Pj)|Xcj) 

+^- 2 (Xcj\Os( 3 S 1 )\x c j) + 0(v 2 T) . (1.2) 

m c 

The term involving the colour- singlet matrix element is the one familiar from the quark- 
potential model. Indeed, up to corrections of order v 2 , one has (Xcj|Ci( 3 Pj)|Xcj) = 
9|Pp(0)| 2 /(27r). The decays of P-wave (and higher orbital-angular-momentum states) 
probe components of the quarkonium wave function that involve dynamical gluons. How- 
ever, the contribution to the inclusive annihilation rate from the higher Fock state |ccg) is 
parametrized through a single number, namely the expectation value of the octet operator 
between the Xcj state, i.e. the colour-octet matrix element in (|1.2j ). 

Consider now the case of exclusive decays of Xcj- For definiteness take Xcj —> and 
define the amplitude Mj through 

T[XcJ - H = — \Mj\ 2 , (1.3) 
m c 

where the aj are real numbers. The amplitude is the sum of a colour- singlet and a colour- 
octet part, Mj = Mj 1 ^ + Mf \ corresponding to the leading and the subleading Fock 
state in (|1.1|), respectively. The singlet amplitude factors into the hard amplitude T^j 
describing the subprocess 

cc 1 ( 3 P J )^qq' 1 ( 1 S ) + q'q 1 ( 1 S ) (1.4) 

and the wave functions for the leading Fock states of the two pions and the Xcj- The Xcj 
wave function can be taken in the non-relativistic limit, in which the charm (anticharm) 
momentum is given by \p ± K and terms linear in the relative momentum K have to 
be kept. The only long-distance information required is hence the singlet decay constant 
/j 1 - 1 or, equivalently, the derivative of the non-relativistic cc wave function at the origin 
in coordinate space, R' P {0) oc fj^y/rn^- The colour-singlet amplitude thus takes on the 
form 
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where if^ is a convolution of the pion distribution amplitudes with a hard kernel, ij = 
J dxdyfa (z)0 w (y) f) 1] (x, y) . 

For the calculation of the colour-octet amplitude two new ingredients enter, the octet 
wave function and the problem of colour conservation. Consider the wave function of the 
|ccsg) state first. It is important to realize that in the |ccsg) Fock state not only the cc 
pair is in a colour-octet state, but also the three particles, c, c and g, are in an S-state. 
Hence orbital angular momenta are not involved and the transverse degrees of freedom 
can be integrated over. Therefore, we only have to operate with a distribution amplitude 
$j(zi, z 2 , z 3 ) (z\ + z 2 + z 3 = 1, where z±, z 2 , z 3 are the plus light-cone fractions of c, c, 



and g, respectively), that is, as usual, subject to the condition J cLzi&^^j = 1, and an 

(8) 

octet decay constant /} for each J. The colour-octet component of the Xcj is then given 
by 

IxS) = f f? J dz.dz^fiz,, z 2} z 3 ) S$ , (1.6) 



where t = A/2 is the Gell-Mann colour matrix and a the colour of the gluon. In the 
following, we will make the plausible assumption that the colour-octet Xcj states differ 
only by their spin wave functions, i.e. the distribution amplitudes as well as the decay 
constants are assumed to be the same for all Xcj states, fj = f^ 8 \ The covariant spin 



wave functions in (|1.6|) are readily constructed 



Sg> = {j> + Mo) (p„/M - 7 „), SS = {i + M a ) . (1.7) 

Owing to the non-relativistic expansion we can employ several simplifications. First, 
the three partons of the |ccg) Fock state can be taken to be collinear up to corrections of 
order v. Then, the c and c three-momenta scale as m c v and differ by at most m c v 2 , i.e. 
Z\ = z 2 up to 0{v 2 ). And finally, the gluon momentum \k\ is peaked at a value of the 
order of the binding energy e = M c — 2m c , where M c is the average charmonium mass. 
Hence we can assume a 5-function-like ccg distribution amplitude, Z\ = z 2 = (1 — z)/2, 
z 3 = z, where z ~ e/M c ~ 0.15. Therefore, analogously to the singlet case, the only 
long-distance information on the octet wave function that enters the final result is the 
octet decay constant f^ 8 \ Identifying the binding energy with m c v 2 (v 2 ~ 0.3), we find 
z 3 ~ v 2 /2 ~ 0.15, in accordance with the above estimate.^] 

Next consider the problem of colour conservation. The obvious solution to colour 
conservation seems to be to demand one of the final-state pions to be also in a higher 
Fock state. The Fock-state expansion of the pion is governed by the ratio of the QCD 
scale and the hard scale Q ~ m c of the process, A = Aqcd/Q- The (multipole-based) 
Fock-state expansion of the pion starts as follows 

|tt(0- + )> = 0(1) IqfcOSo)) + 0(A) \qq & { 1 Px)g) 

+0(A 2 ) {|qq 8 ( 3 5 1 )g) + ...} + 0(A 3 ) . (1.8) 

Both octet terms give corrections of the same order A 2 with respect to the leading one: 
the extra suppression of the magnetic dipole transition for 3 S\ (A 2 rather than A as for an 

^ince we are working to lowest order in the w-expansion we may simply take p c .c = Zi^p, k = z^p in 
the actual calculation; the amount of off-shcllness is p\ — ~ m^v 2 and k 2 ~ m 2 v A . Correspondingly 



the gluon has three polarization states, cf. (1.7) 
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electric dipole transition) is compensated by the suppression oc A of P-wave production. 
The octet amplitude will then look as follows: 



f f (p,8) m Am 

M(8)rl T3 1 2 1 \ Jt J n J 1 J 
y[ P x \ ~m c a s m c - 

m c m A c m l c m c \ z 
f f (5 ' 8) f (8) /i S ' 8) 

M(8)r3ol 2/ \ Jir Ji J J J 



1.9) 



where /} is a convolution of the hard process 



cc 8 ( 3 ^) ^qq^( 1 P 1 ,^ 1 )+q'q 1 ( 1 S' ) , (1.10) 

the ordinary pion distribution amplitude and the distribution amplitude for the higher 
Fock states of the pion. We have indicated the extra suppression of these higher Fock 
states by assuming the octet decay constant to scale as f^ p ^ ~ \f w , fi s ' 8 ^ ~ A 2 /^, and 
indicated explicitly the extra suppression of P-wave production. The factor A 2 (oc 1/m 2 ) 
appears as a consequence of this particular solution of colour conservation: the constituent 
gluon merely acts as a spectator which runs from the Xcj to one of the pions without 
changing its momentum. 

Yet, we disregard this possibility of colour conservation for the following two reasons. 
First, it requires the specification of two new distribution amplitudes for the two higher 
Fock states of the pion. Second, the Fock-state expansion (|1.8| ) for the pion might be 
badly convergent. Even if a Fock-state expansion existed, it need not obey the usual 
multipole expansion assumed in ( |1.S| ). If we do not want to work with explicit higher Fock 
components of the pion, we have to answer two questions: what the constituent gluon of 
the |ccg) state couples to, and what determines the probability, i.e. what the analogue of 
the charmonium /( 8 ) is on the pion side. To this end we invoke two assumptions. 

First, we take QCD perturbation theory to be valid down to virtualities of the order 
of zm 2 ~ (m c v ) 2 (see Figs. 4 and 5 and Appendix). Then we can attach the gluon of the 
|cc 8 g) state (in all possible ways and with a coupling a s ° ft ) to the hard process leading to 
the Feynman diagrams shown in Figs. (| and [|. Thus the hard process is 

cc 8 ( 3 ^)g - qqU 1 ^) + q'q^So) • (1-11) 

Second, the transverse motion of the valence quark and antiquark determines the 
importance of the momentum distribution of the constituent gluon. It is clear that the 
diagrams corresponding to ( |1 . 1 1| ) must contain contributions, which, to leading order in 



a s and z, constitute the two higher Fock states |qqg) of the pion. In the approximation of 
collinear light quarks, one thus encounters singular integrals. These infrared singularities 
precisely correspond to the long-distance wave functions describing the higher Fock states 
of the pion. In a collinear calculation one hence needs distribution amplitudes for the 
higher Fock states into which the singularities can be absorbed. In our ansatz we keep 
the (non-perturbative) transverse motion of the pion's valence constituents (quark and 
antiquark). Then all integrals are finite. More precisely, gluon momentum configurations 
are selected such that the q and q relative momentum "matches" the one prescribed 
by the pion wave function. In this sense one can say that our approach leads to the 
dynamical generation of higher Fock components. We emphasize that also in this scheme 
of colour conservation the colour-octet contribution is not suppressed by powers of l/m c 
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as compared to the colour-singlet one. This has been discussed by us in [|J, see also 
Sect. 4. 

In our previous calculation ||, we employed the collinear approximation and simply 
regularized the singular integrals through a cut-off, which corresponds to an average 
transverse momentum of the quarks inside the pions. The final results for the colour- 
octet contribution to the decay amplitude then basically depended on a single parameter 

K = yafft/^/f? 2 • (1.12) 

Here af {t denotes the coupling of the gluon of the |ccg) Fock state to the hard process, 
and is the octet decay constant. There is an additional contribution from terms not 
singular for g — > 0, which is proportional to \J(yl° ft f^- It is, however, rather small, less 
than about 10%. 

In this work we are going to extend our previous one by working within the mHSA. 
This leads to three major improvements. First, the scale of the coupling constant is not 
fixed in the collinear approximation, implying a big uncertainty. In contrast, the scales 
of a s are fixed in the mHSA. Second, in the collinear approximation there is usually no 
suppression of the end-point regions of the distribution amplitude, while these infrared 
regions are explicitly suppressed in the mHSA through a Sudakov form factor. And 
finally, no ad-hoc regularization of singular integrals appearing in the calculation of the 
octet contribution is necessary. Rather, the momentum distribution of the octet gluon is 
linked to the k± dependence of the pion wave function. 

The paper is organized as follows: In Sect. 2 we present our ansatz for the pion wave 
function and in Sect. 3 we introduce the mHSA in the case of the colour-singlet decay 
contribution to XcJ 7r + 7r~. In Sect. 4 and Appendix, a detailed description of the 
colour-octet contribution and a discussion of the results for the decay widths are given. 
In Sect. 5 we present results for the Xcj decays into pairs of kaons and etas. Finally, in 
Sect. 6 we comment on J/ip decays, before we end our paper with a summary and some 
conclusions (Sect. 7). 



2 The pion wave function 

We write the pion's |qq ) Fock state, which is the leading one, in a covariant fashion 

where the spin wave function is defined as (the pion mass is neglected throughout) 

S, = A=H ■ (2-2) 

Integrating the light-cone wave function ty n over k±, the transverse momentum of the 
quark with respect to the pion momentum, up to a factorization scale /ip, one arrives, 
up to a constant dimensionful factor f 7T /(2 v6), at the distribution amplitude ^(x, fxp). 
The constant factor plays the role of the configuration space wave function at the origin; 
f n = 130.7 MeV is the usual pion decay constant. 
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Upon expansion over Gegenbauer polynomials C^/ 2 \ the distribution amplitude at a 
scale \ip can be characterized by non-perturbative coefficients B n (see [[!]] and references 
therein) : 



1+ J2 Bn^0) 



n=2,4,. 



' a s (ld F y 

«s(/io), 



(2.3) 



The odd expansion terms do not appear since, in the isotopic limit, the pion distribution 
amplitude is symmetric (f) w (x) = </> n (l — x). In ( |2.3|) a s is the strong coupling constant, 
and fx a typical hadronic scale, 0.5 1 GeV. Since the anomalous dimensions 7 n in 

Q2.3|) are positive fractional numbers increasing with n, higher-order terms are gradually 
suppressed and any distribution amplitude evolves into 0as(^) = Gx(l—x) asymptotically, 
i.e. for ln(Q 2 /AQ CD ) — > oo. The asymptotic (AS) distribution amplitude itself shows no 
evolution. 

From the investigation of the pion-photon transition form factor F yn (Q 2 ) ||, ^ || |H], 
lT| it follows that the form of the pion distribution amplitude is very close to the asymp- 
totic form. Also recent QCD sum-rule analyses |12]| lead to that result. Hence, all terms in 
the expansion ( |2.3|) of the pion distribution amplitude with n > 2 will provide only small 
corrections to exclusive observables and it is legitimate to truncate that expansion at the 
second term (note that the frequently used Chernyak-Zhitnitsky distribution amplitude 
p~3|1 is given by B 2 = 2/3 and B n = for n > 4) and consider now B 2 as the only soft 
parameter. 

In the mHSA we consider the Fourier transform of ^ n to transverse coordinate space 
which, following ||, [7J, is written as 



2VQ 



(2.4) 



The dependence on the transverse separation b, canonically conjugated to k^, is thereby 
chosen to be of a simple Gaussian form 



b, fi F ) = "In exp 



x(l — x) b 2 
4a2(/i F ) 



(2.5) 



Here, the momentum fraction x and b (or kjj refer to the quark; the antiquark momentum 
is characterized by 1 — x and b (or kj_) throughout. The pion's transverse size parameter 



<V = cltt^f) is fixed from the process n° — > 77 [fP|] . That constraint leads to the 
closed formula 1/a 2 = 8 (1 + B 2 {^f)) tt 2 fl under the assumption B n = for n > 4 
(for the asymptotic distribution amplitude (B 2 = 0) a n = 0.861 GeV -1 |4]]). The scale 
dependence of a n is an approximation sufficient for our purpose. Using ( |2. 1| ) ( p~5|) we 
obtain the following expressions for the probability of finding the pion in its valence Fock 
state, for the mean transverse momentum and the mean radius of the qq Fock state 



P - — - 

^qq 4 



(kl) 



1 1 + 18/7-Bj 
1 + B 2 

fl (1 + B 2 



R, 



qq 



47T 2 

5 

3 
8^2 



1-6/7^2 + 12/7 B, 



1 + 18/7SI 



1 + 3B 2 + 5A/7B 2 
(l + B 2 y 



(2.6) 
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Figure 1: B 2 dependence of the probability P qq - (upper left), the radius R (upper right), 
r.m.s. transverse momentum (k\) l l 2 (lower left) of the pion's valence Fock state as well as 
the shape of the distribution amplitude (lower right). In the lower right figure, the solid 
(dashed) line corresponds to B 2 = (—0.2). 



As shown in Fig. our wave function has appealing features as a function of B 2 . For 
the asymptotic form of the distribution amplitude we obtain the values P qq - = 0.25, 
(/c^) 1 / 2 = 367 MeV and P qq - = 0.42 fm. We remark that the latter is considerably smaller 
than the pion's charge radius of 0.66 fm to which all the Fock states contribute. Since 
P q q is much smaller than unity, higher Fock states are important components of the pion 
and, of course, have been seen for instance in Drell-Yan dilepton (/i + /x~) production in 



p ion-proton collisions fL6fl . At large Bjorken x, where the valence quarks dominate, our 
wave function, for small B 2 , is consistent with the valence quark distribution function 
as derived in [M (see [|]]). For large negative values of B 2 , i.e. at a very low scale, the 
wave function bears resemblance to a constituent wave function with a probability close 
to unity and a valence Fock state radius that is approximately equal to the charge radius 
of the pion. 

In exclusive reactions, contributions from higher Fock states of the pion are usually 
neglected since they are suppressed by inverse powers of the relevant hard scale, see ( |1.8|) . 
As explained in the introduction, we do not make use of explicit Fock-state wave functions 
of the pion. Rather we generate the contribution corresponding to a constituent gluon of 
the pion dynamically. This is achieved in two steps: (i) the constituent gluon of the Xcj is 
coupled perturbatively (oc a s s oit ) to the hard process; (ii) the dependence of the momenta 
of the pion's valence constituents (i.e. the q, q momenta) on k± is kept. Since it is the 
non-zero value of k± that regularizes the propagators associated with the "would-be" 
constituent gluon of the pion, the k± dependence of the wave function of the qq valence 
Fock state determines the size of the qqg contribution in the pion. 
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3 The colour singlet contribution and the mHSA 

In [|[] we have analysed the colour singlet contribution to Xcj — ► tt + 7i~ already within the 
mHSA. Here, we will give a more detailed description of the calculation. 

The colour singlet components of the Xcj are expressed in terms of non-relativistic 
wave functions 

8a f d 3 k 



Y (l)\ _ f dfc jf,(l) f M c(i) 



v/3 

where the represent reduced wave functions, i.e. full wave functions with a factor of 
the relative momentum removed from them (Kp = and, for instance, = (0, k) 
in the meson rest frame). The denote the covariant spin wave functions 

SjP = -L + M + 2Jfl If, 



c(l) _ 1 

2 "71 



2 



+ m 2 ) 7p + — + m 2 )k p -m % 



2 



e»°K a . (3.2) 



The spin wave functions represent an expansion upon powers of K up to terms of 0(K 2 ) 
|T7| The derivative of the non-relativistic cc wave function at the origin is introduced 
by 



R' P (0) = i ^ ^ ^T c [ ^rr^ = ^Vl&;/o (1) - (3.3) 
3 J lo7r s Mo 

In the non-relativistic approximation one has Mq ~ M 2 — 2m c and the same normal- 
ization of the two wave functions (usually unity) Pq c = P 2 C . In this case one finds 
^2 = V^^o an d the same relation between the two singlet decay constants fj. 

The starting point of the calculation of the colour singlet decay amplitude within 
the modified HSA is the convolution with respect to the momentum fractions x, y and 
transverse separation scales b 1; b 2 of the two pions 

„, m/ , , 32 v^tt 372 |^p(0)| f 1 , , f d 2 b x d 2 b 2 



3V3ml /2 Jo " J (4tt) 2 (4tt) 2 

%(y,b 2 )f HJ (x,y,b 1 ,b 2 )^ 7T (x,b 1 ) exp[-S(x, y, b u b 2 , t x , t 2 )\ , (3.4) 

which adapts the method proposed by to our case of exclusive charmonium decays 
((T = l,cr 2 = v3/2). T HJ is the Fourier transform of the hard-scattering amplitude T HJ 
with the fc^-dependence retained 



f HJ (x,y,b 1 ,b 2 ) = f \\ \l T HJ (k ±1 ,k ±2 ) exp [-«kj. r bi - 2k ±2 -b 2 ] . 



T HJ is to be calculated from the graphs shown in Fig. |2] and reads 



2 The spin wave functions can be written in a more compact form: Sj = Sj + {S^, If} /Mj where 
j 1 ^ represents the 0(K) spin wave function. 
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Figure 2: Feynman graphs for the colour- singlet decay XcJ — ► ^-n (J = 0, 2). 



T HJ (x,y, k ±1 ,k ±2 ) 



Oi B (ti) a s (t 2 ) 



1 + 



(-2) J / 2 (x-yf 



(g 2 + ie)(g 2 2 + ie)N \~ ' 2 
With 

K ± = kj.i-kj.2 
the virtualities of the internal quark and the two gluons are 

N = x(l - y) + (1 - x)y + K 2 ± /(2m 2 c ) , 



N 



(3.6) 



(3.7) 



~g\ = xy- K 2 J(4m 2 c ) , 
~g 2 2 = (l-x)(l-y)-K 2 J(4m, 



2 ) 

cJ ' 



(3.8) 



where a common factor 4m\ 1S pulled out. The convolution formula (3.4) involves two 
independent transverse separation scales b\ and b 2 , each associated with one of the two 
pions. The foj also provide the factorization scales hfi = (i = 1,2) of the pion 
wave functions. Distinct from these scales are the renormalization scales tj (j = 1,2) 
corresponding to the virtualities of the hard gluons. 

The fact that T HJ depends on k ± i and kj_ 2 only in the combination implies the 
following result for the Fourier transform of the hard-scattering amplitude {b = 2m c bi) 



T HJ {x,yMM) =Am 2 c a s {tl)a s {t 2 2 )5^ 2 \h 1 - b 2 ) 



-i/4 



H 



xy b) 



x + y 



1 + 



(1 -x)(l -y)-xy 
(-2y/ 2 (x-y) 2 \ 



2(x + y) 



H 



l-x)(l-y)b) 



x-y 



1 + 



-2)H\x - yf 
2{2-x- y) 



+ 



lK U(x(l-y) + (l-x)y)/2b) 



1 + 



7T 



+ 



8vr 



(-2y/ 2 (x-y) 2 
(x + y)(2-x-y) ^ (x + y){2 - x - y) j 

(_2)^(x - yf bK^jxjl - y) + (1 - x)y)/2b) \ 

(x + y)(2-x-y) ^J( x (l-y) + (l-x)y)/2 J 



(3.9) 
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Hq 1 ^ and Kj denote Hankel and modified Bessel functions, respectively. In physical terms 
the appearance of the ^-function, which simplifies the numerical work enormously, means 
that the two pions emerge from the decay with identical transverse separations. 

The novel ingredient of the mHSA is the Sudakov factor exp[— S], which takes into 
account those gluonic radiative corrections not accounted for in the QCD evolution of the 
wave function. In next-to-leading-log approximation the Sudakov exponent reads 



S(x,y, bi,b 2 ,ti,t 2 ) = s(x,bi,2m c ) + s(l - x, h,2m c ) 

+ s(y, K 2m c ) + S{ l-y, K 2m c ) - 1 log , T^P^M^ 1 » 

P log(l/(6lA Q CD)) log(l/(& 2 AQCD)) 



(3.10) 



where the function s(x,b,Q), originally derived by Botts and Sterman || and later on 
slightly improved, can be found, for instance, in The last term in (3.10) arises from 
a renormalization group transformation from the factorization scales to the renormal- 
ization scales tj at which the hard amplitude Th is evaluated. The renormalization scales 
tj, entering Tjj as the arguments of the two strong coupling constants a s , are chosen as 



ti = max{4xym 2 c , 1/bj, 1/bQ , t 2 = max{4(l - x)(l - y) m 2 c , l/b\, l/b\} 



(3.11) 



thus avoiding large logs from higher-order pQCD. In the limit b\ — > 1/Aqcd the tj may 
in principle approach Aqcd at which the one-loop expression for a s is diverging. Still 
the integral appearing in (3.4) is regular since the Sudakov factor compensates the a s 
singularities. Actually the suppression of the end-point regions is so strong that the 
bulk of the perturbative contribution is accumulated in regions of small values of a s . 
The physical picture behind the Sudakov suppression is that qq pairs with large mutual 
separation tend to radiate so many gluons that it becomes impossible for the hadronic 
state to remain intact and that the exclusive process can take place. 

Replacing exp[— S] by 1 and ignoring the transverse momenta in Tjj, one finds from 
(3.4) the colour singlet decay amplitude within the standard HSA as derived by Duncan 
and Mueller |19| . In that approach the renormalization scale is taken as the charm quark 
mass and customarily identified with the factorization scale. 

In terms of the amplitude (3.4) the Xcj decay widths into pions are given by 



r[xco 
r[xc 2 



7T 7T 



7T 7T 



32 irm c 
1 



M( X co - 

M{ X c2 



7T 7T 



240vrm„ 

The results obtained within the mHSA can be cast into the form 



J 7T 



m° 



r[Xcj - Tr+Tr-] = ^ |i^(0)| 2 « s (m c ) 4 a { p + bf BM + c [ j } B 2 (fM 



(3.12) 



(3.13) 



where the coefficients aj, bj, and cj are complex- valued. Written in the form ( [3.13] ) we 
may then immediately find the prediction for a distribution amplitude with a given value 
of the expansion coefficient B 2 at the scale /xq- It is still convenient to divide out the fourth 
power of a s at the fixed scale m c in Q3.13 ), since the main effect of the strong coupling is 
thus made explicit. The actual effective renormalization scale /ir (i.e. the mean value of 
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(Ir) in the modified HSA differs from m c and depends on B 2 (e.g. = 1.15 GeV and 
hence a s = 0.43 for -82(^0) = 0). 

In the following we will choose as central value -Rp(O) = 0.22 GeV 5 / 2 and m c = 
1.5 GeV, which is consistent with a global fit of charmonium parameters |20| as well as 



results for charmonium radii from potential models [21]. Moreover, we use the one-loop 



expression for a s with Aqcd = 200 MeV and four light flavours. As we have demonstrated 
in H by varying the input parameters m c and R' P (0), B 2 , and Aqcd, the colour-singlet 
contribution alone, i.e. a calculation based on the assumption that the Xcj is a pure cc 
state, is insufficient to explain the observed decay rates. 



4 The colour octet contribution 

Recently, the importance of higher Fock states in understanding the production and the 
inclusive decays of charmonia has been pointed out ||. The heavy-quark mass allows 
for a systematic expansion of both the quarkonium state and the hard, short- distance 
process and, hence, of the inclusive decay rate or the production cross section. The 
expansion parameter is provided by the velocity v of the heavy quark inside the meson. 
The crucial observation is that, for inclusive decays (and also production rates) of the 
Xcj-, both states in fll.ip contribute at the same order in v and, hence, the inclusion of the 
"octet mechanism", i.e. the contribution from the |ccsg) state, is necessary for a consistent 
description. (Without its inclusion the factorization of the decay width into long- and 
short-distance factors is spoiled by the presence of infrared-sensitive logarithms.) The 
inclusion of the octet mechanism is also necessary in exclusive charmonium decays in 
order to get a consistent description. This is so because the usual suppression of higher 
Fock states in exclusive reactions does not appear for the decay of P-wave charmonium: 
it is compensated for by the P-wave nature of the cci Fock state M. 

As a consequence of employing the collinear approximation we encountered in || 
a number of singular integrals, which we regularized by a parameter g related to the 
(neglected) transverse momentum of the internal quarks and gluons. However, within 
the modified HSA where the dependence on the partonic transverse momenta is taken 
into account, all integrals are finite and there is no longer any need for the parameter 
g. Apart from that, the modified HSA also provides an explicit prescription for the 
renormalization scales entering the strong coupling constant such that a s is no longer a 
free parameter. Altogether, the modified HSA not only guarantees a self-consistent and 
thus reliable calculation of the perturbative contribution, it also helps in the present case 
to diminish uncertainties of the collinear approximation. 

As explained in the introduction, the Fock-state gluon has a typical momentum frac- 
tion z ~ m c v 2 /2m c of the Xcj momentum. We treat z as a small number. We have 
checked that the results are nearly ^-independent for 0.1 < z < 0.3. Hence the |ccg) 
wave function <&f of Q can be taken as exhibiting a 5-function-like peak at the value 



z 3 = z. Moreover, in our results we will neglect all terms of order z 3 and higher. This 
is in line with the non-relativistic expansion of the charmonium decay rate since z scales 
as v 2 . Apart from z, the only free parameter on the charmonium side is f( 8 \ the colour 
octet wave function at the origin. 

The details of the calculation of the colour octet contribution within the modified HSA 
can be found in the Appendix. We have divided the diagrams contributing to the colour 
octet decay into eleven groups, numbered by the index i. The contribution of each group 
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is of the form (3.4). Furthermore, we have introduced renormalization scales tij, which 
are determined by the virtualities of the internal quark and gluon lines (see Appendix) 
and depend non-trivially on the integration variables. As already discussed in Sect. [5] 
we may simply evaluate the running coupling constant by its one-loop expression at the 
scales because its singularity (to be reached for b\ or b 2 - ► Aqcd) is compensated by 
the Sudakov factor. 

We now want to direct the reader's attention to the question of gauge invariance. 
Treating the problems as presented above, the results we obtain are gauge invariant to 
order z 2 (with z being the fraction of the Xcj momentum carried by the valence gluon), 
which is sufficient, because we neglect terms of 0(z 3 ) in our results. Eventual violations 
of gauge-invariance at order z 3 may be traced back to the general problem to describe 
a constituent within the parton model. This conjecture is supported by the following 
observation we made: if one changes the c-quark propagator mass in groups 8-10 from 
the usual value m c to (1 — z)m c , we end up with completely gauge-invariant results. 
Admittedly, this choice is by no means mandatory and we do not apply it actually, but 
it reveals that the violations of gauge invariance to 0(z 3 ) are related to the neglect of 
binding effects for which that particular choice of the mass in the c-quark propagators 
seems to compensate for. 

In order to combine colour singlet and octet contribution we expand the colour-octet 
decay amplitude (A.l) in a fashion analogous to the colour singlet case: 



4 



y/XjTrm c I R'JO) I \a s (m c ) fir 



M^( XcJ - vr+O = (af + bfB 2 + of Bl) 



(4.1) 



where Ao = 32 and A2 = 240. Thus, to obtain the sum of singlet and octet contributions, 
we simply have to add the coefficients, e.g. aj 1+8 ' ) = a^p + af\ 

First, we determine a value for the unknown decay constant f^ 8 \ which is contained 
in the coefficients , bj and Cj' ) , from a fit of our predictions to the experimental data 
on the decay widths, using the asymptotic pion wave function (B 2 = 0). We obtain 

/ (8) = 1.46 x 10" 3 GeV 2 , (4.2) 

which is larger than the result estimated in Ref. |§. This is mainly a consequence of 
the treatment of the strong coupling constant a s . In Ref. || we have chosen 0.45 for 
the coupling of the two hard gluons and assumed al oit ~ ir for the coupling of the soft 

Xcj Fock-state gluon to the hard process: a s (ti)a s (t 2 ) a s (t3) = 0A5 2 ^ ~ 0.36. In 
the present analysis based on the mHSA we may extract effective values for the strong 
coupling constant. We distinguish the hard a B , evaluated at the scales ta and from the 
soft one at t^. Using the asymptotic form of the pion distribution amplitude we obtain 



. ft , ,„ 0.492 2 x V0^7 = 0.181 in X o -> 

(a s tod ) («s ) 7 = I X . (4.3) 

[ 0.509 2 x t/oIqE = 0.195 in X 2 -> tt+tt- 

The smaller values of the products of strong couplings have to be compensated for by 
a larger value of f^ 8 \ Observe that the coupling o^ oft of the constituent gluon of the 
Xcj is only slightly larger than the coupling of the "hard" gluons. This is caused by two 
effects: (i) the dynamical setting of the scales makes a^ ard larger than the naive, collinear 
estimate a£ ard = a!? ard (m c ); (ii) the Sudakov factor is effective enough to suppress the 
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J = 


J = 2 


CI \ 

af 


23.36 + 14.67? 


5.14 + 3.42 1 


(8) 

a y 


33.90 + 15.89 z 


11.60 + 4.12? 


a J 


57.26 + 30.56 % 


16.74 + 7.542 



Table 1: Contributions to aj (singlet, octet, and their sum for /( 8 ) = 1.46 x 10 3 GeV 2 ). 

soft phase-space regions, in which af ft would become large. The fact that this coupling 
is not probed at large values gives us confidence in our perturbative treatment of the XcJ 
cosntituent gluon. 

The question remains whether the result (|4.2j ) is physically sensible or not. Therefore, 
we are going to estimate the probability P C c g to find the Xcj in the colour-octet state from 
the following examplary parametrization for the ccg wave function 

y { o\zi, k±i) = Nz x z 2 zl exp {-2a\m\ [{z 3 - zf + {z x - z 2 ) 2 ] } 

xexp{-a 2 x J2 k li} , ( 4 - 4 ) 

This ansatz combines the known asymptotic behaviour of a distribution amplitude for 
a qqg Fock state with a mass-dependent exponential and a Gaussian k± dependence in 
analogy to the Bauer-Stech-Wirbel parametrization of charmed-meson wave functions 
p2fl . The mass- dependent exponential guarantees a pronounced peak of the distribution 
amplitude at Z\ ~ z 2 — (1 — z)/2. The 5-function-like distribution amplitude used in the 
estimate of the colour-octet contribution appears as the peaking approximation to this 
function. Since the ccg Fock state is an S'-wave state we assume its radius to be equal 
to that of the S'-state charmonia, namely 0.42 fm ||21|| . In this case the transverse size 
parameter a x takes the value 1.23 GeV -1 . The probability of the colour-octet Fock state 
is then found to be 

P c - Cg = (/ (8) /2.1 x 10" 3 GeV 2 ) 2 . (4.5) 

This relation is obtained with z = 0.15; it is however practically independent of the exact 
value of z. We stress again that the ansatz ( |4.4|) is only an example for a possible wave 
function. Thus, our result of about 1/2 for P C c g (from ( |4.2| ) in ( |4.5| )) is not in apparent 
disagreement with expectation and we may conclude that the result ( |4.2|) corresponds to 
a large but not unphysical value for P C c g . 

Using the value fl4.2| ) the coefficients aj add up as shown in Table |l| Our results for 
the decay widths are presented in Table ||. Note that in the colour octet decay amplitude 
there are terms contributing to decays into neutral pions only (see group 11 in Fig. |5|). 
We may thus find deviations from 1/2 for the ratio of the decay width into neutral pions 
over that into charged pions. In experiment as well as in our calculation these deviations 
are found to be small. As discussed before we perform a peaking approximation in the 
momentum fraction z carried by the Xcj Fock-state gluon. As in Ref. |8| we set z = 0.15. 
We have checked that our results are only weakly dependent on the actual value of z 
within the range between 0.1 and 0.2 and all our conclusions remain valid. 

In Fig. |3] we sketch the dependence of the Xco ~~ > 7r + 7r~ decay width on the expansion 
parameter B 2 in the allowed B 2 range [Hj. We dispense with a corresponding plot for x c2 
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PDG [23| 


BES [24] 


r[Xco - vr+TT-] [keV] 


45.4 


105 ± 47 


64 ±21 


T[Xc2 -> 7T + 7T-] [keV] 


3.64 


3.8 ±2.0 


3.04 ±0.73 


r[Xco^vr°7T ] [keV] 


23.5 


43 ± 18 




r[Xc2^7r°7r°] [keV] 


1.93 


2.2 ±0.6 





Table 2: Results for the XcJ decay widths into pions (/ (8) = 1.46 x 10~ 3 GeV 2 ; B 2 = 0) 
in comparison with experimental data. The BES result for T[x c o 7r + 7r~] is evaluated 
with the BES result for the total width. In the other cases the PDG average values for 
the total widths are used. 




Figure 3: Dependence of the prediction for the Xco ~~ > 7r+7r decay width on the expansion 
parameter B 2 of the pion distribution amplitude. 
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since its width is found to depend on B 2 in the same way as the width of Xco ~^ k+h^ . 
As can be seen from Fig. § the dependence on B 2 is quite strong, i.e. pushing B 2 from 
the central value zero to its upper limit B 2 = 0.1 doubles the width. If the octet decay 
constant is determined by some other method, the XcJ — > decays constitute a severe 
test of the pion distribution amplitude. 

We conclude from these considerations that it is possible to explain the experimental 
data on XcJ — ► 7r + 7r~ within a perturbative approach, using the asymptotic form of 
the pion distribution amplitude, provided colour-octet contributions are included. The 
numerical results of our new calculation within the mHSA are similar to our old ones 
||, obtained within the collinear approximation. The new calculation within the mHSA 
improves, however, the old one in various respects: (i) the scales of the coupling constant 
are determined, (ii) soft phase-space regions are suppressed, and (iii) all propagators 
associated with the constituent gluon are regularized through the non-zero transverse 
momenta of the pion's constituents. 



5 XcJ decays into other light pseudoscalars 

The considerations made so far for the case of the pion can straightforwardly be gener- 
alized to the decays of Xcj into the other light pseudoscalar mesons, e.g. etas and kaons. 
The r\ meson is a linear combination of an SU(3) octet and singlet component 

| r] ) — cos dp | r/ 8 ) — sin 9 P | rji } . (5.1) 



For the wave functions of rji and rjs we follow [|J] and make the same ansatz (|2.1|) - (|2.5|) for 
it as for the pion (with a v = a^). The decay constants f vi and f m as well as the mixing 
angle 6p have been determined from a fit to the 77-7 and rj'-^f transition form-factor data 



f m = 145 ± 3 MeV; f m = 136 ± 10 MeV; 9 P = -18° ± 2° . 

(5.2) 

Since the strong interaction is flavour-blind, the hard scattering amplitudes for the decays 
of Xcj into rji and rjs pairs equal both that of the 7r 7r° channel (with group 11 included), 
i.e. the only modification with respect to the pion case is to replace / 2 by cos 2 #p/ 2 g + 
sin 2 #p/ 2 = (144 MeV) 2 in Eq. (|3.13|) . In addition we take into account the r\ mass in 



the phase space factor, [1 — m 2 //??. 2 ] 1 / 2 , of the decay width (|3.12|) , ( |3.13| ). Using, as in 



B 2 = for the rj distribution amplitude, we obtain the Xcj decay widths into pairs of etas 
as quoted in Table |3|. With that asymptotic form of the rj distribution amplitude, fair 
agreement with the data is achieved. On the other hand, a fit of B 71 to the T[x c j w] 
data yields a value of -0.036 for it and leads to a slightly better agreement with the data. 

The kaon state is also written in the form (|2.1| ) - fl2.5|) , with /„- exchanged for fx = 
157.8 MeV. Again in order to simplify matters, we assume ax = cl-k- Since the kaon 
consists of up and strange quarks with different masses one may no longer expect the 
kaon distribution amplitude <pK to be symmetric under exchange x <-» 1 — x. Therefore, 
one has to include in the expansion Q2.3|) also the antisymmetric terms 



i+ t *<*)(*grc?«(2,-i: 



71=1,2,. 



(5.3) 
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r[ XcJ ^/i+K-] [keV] 


T[Xcj tjij] [keV] 




J = 


J =2 


J = 


J = 2 


flf = -0.176 


22.4 


1.68 






b? = -0.100 


38.6 


2.89 






B\ = -0.036 






24.0 


1.91 


Bl = 0.000 






32.7 


2.66 


PDG @ 


99 ±49 


3.0 ±2.2 


35 ±20 


1.6 ± 1.0 


BES [24] 


52 ±17 


1.04 ±0.43 







Table 3: Results for the Xcj decay widths into kaons and etas in comparison with exper- 
imental data (/W = 1.46 x 10~ 3 GeV 2 , Bf = 0). 

where e.g. the anomalous dimension 71 is 32/81 [|TJ. 

Information on the antisymmetric part of the kaon distribution amplitude may be 
extracted from the valence quark distribution functions of the kaon at large x, where the 
contributions from Fock components higher than the valence Fock state are negligible. 
In particular the ratio of strange antiquarks over up quarks is sensitive to the amount of 
asymmetry in the kaon distribution amplitude. Unlike the cases of the nucleon and the 
pion, no systematic model-independent analysis of the kaon quark distribution functions 
exists so far. From calculations of the kaon's valence quark distributions within the 
Nambu-Jona-Lasinio model [25f|, there is evidence that at x ~ 0.8 there are about twice 



as many strange antiquarks as u quarks inside the K + . Truncating (|5.3| ) at n = 2 one 
sees that such a ratio can be accommodated by a coefficient Bf of order —0.1. About 
the same value of By is found in an instanton model p6| . 



Due to the symmetry of the hard scattering amplitude for Xcj ~ > K + K under the 
simultaneous exchanges x «-> 1 — x and y <-> 1 — y, B± will enter the decay amplitude 
only quadratically, i.e. in ( |3.13|) and ( |4.1|) the only additional terms are d^j\B^) 2 . For 



[By I < 0.1 the contributions from the antisymmetric part of the distribution amplitude 
therefore provide only tiny corrections to the decay of Xcj into kaons and can be neglectedP]. 
Hence, the hard scattering amplitude for Xcj decays into K + K~ has the same form as 
the one for the decays into charged pions. A inspection of Table [| reveals that already a 
small negative B^ value of the order of —0.1 to —0.2 is sufficient to obtain fair agreement 
between the phase-space-corrected mHSA result and the available data. Admittedly, a 
precise value of B^ cannot be determined at present because the BES data and the 
PDG average values only agree within large errors; we therefore present results for two 
different B^ values in Table |3|. Despite this uncertainty it seems that the kaon distribution 
amplitude is somewhat narrower than the asymptotic one (which is favoured for the pion) . 
This finding appears to be plausible considering the comparatively large strange quark 
massfU For B± = 0, B^ = —0.1 we obtain from ( |2.6| ) the following properties of the 

3 This was already observed by the authors of 
Remember that the non-relativistic case of infinitely heavy quarks is described by a (5-function-like 
distribution amplitude. 
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kaon's valence Fock state: probability 0.37, r.m.s. transverse momentum 381 MeV, and 
radius 0.44 fm. 

6 Comments on J/ip —> tt + tt~ decays 

It has long been known that there exists a problem with the pion form factor. The value 
sF n = 0.95 ± 0.08 GeV 2 at s = M)^ deduced from the decay J/if) — > 7r + 7r is much larger 
than the space-like one at the equivalent Q 2 scale (Q 2 = —s): Q 2 F 7T = 0.35 ± 0.10 GeV 2 . 
Admittedly, the space-like data |27j suffer from large systematical errors. The discrepancy 
is even more severe for ip(2S) — > tttt (sF n = 2.67 ± 0.87 GeV 2 ). Unfortunately, the data 
on the time-like form factor measured in e + e~ — > 7r + 7r~ [Q, suffer from low statistics and 
are inconclusive. While, within the large experimental errors, these data are in agreemnt 
with the value for the time-like pion form factor as extracted from J/ip —>■ tttt, a much 
smaller value (close to that of the space-like form factor) cannot be excluded. These two 
differences between the time-like and space-like pion form factors pose a true challenge 
for the HSA. Within the standard HSA there is no difference between the form factor in 
the two regions (except for interchanging Q 2 with — s). While possible explanations for 
a larger value of the form factor in the time-like region than in the space-like one within 
the mHSA have come up ||29|| , the large values of the pion form factor extracted from 
ip{nS) — > 7T7T still lack a decent explanation. Here we want to propose a solution to this 
puzzle. 

The one-to-one correspondence between sF^ and J/ip — > ttti is based on the conven- 
tional assumption that in a leading-twist analysis the only contribution to J/ip — > mi 
comes from the electromagnetic process with one intermediate photon. The hard process 
is 0{a 2 a 2 m ). Purely hadronic processes, seemingly larger by factors ct p s /a 2 m) are usually 
supposed to be suppressed by powers of Aq CD /m c . Moreover, the rate for the 0{aP s ) 
contribution via three intermediate gluons is zero if the light-quark masses are assumed 
to be zero f3"0fl . 

Recent developments [[J in the theory of heavy quarkonia have shown that the appro- 
priate expansion parameter is not Aqcd/^c, as is the case for light systems, but rather 
the velocity v of the heavy quark (or antiquark) in the bound state. The dominant Fock 
state of the J/ip is |(cc)i( 3 5'i)), a colour-singlet cc pair in a spin-triplet S-wave state. This 
state decays to a pair of pions via a virtual photon and its rate is of order a 2 m a 2 v 3 . (The 
factor v 3 arises from the squared wave function). Contributions from higher Fock states 
are suppressed by powers of v and af . Factors of v give the probability to find these 
soft gluons in the quarkonium, while factors of a s s oit are associated with the coupling of 
the soft gluons to the decay process. 

Corrections to the electromagnetic decay of the valence Fock state first arise at relative 
order v 4 from the |(cc) 8 ( 3 S'i)gg) Fock state and the |(cc) 8 ( 3 -Pj)g) Fock state. In both cases 
the hard process is of order oq, while the soft part scales as v 3 (v 2 a s s oit ) 2 (v 3 v 2 (v 2 a s s oit )) 
for the former (latter) Fock state. (The former contains two soft gluons, the latter only 
one, but it is a P-wave annihilation. For comparison, in the conventional HSA (all gluons 
are purely perturbative and the valence Fock state cc"i( 3 Si) is the only relevant one), the 
diagrams corresponding to |(cc) 8 ( 3 5'i)gg) are of order a®. The rate for the 0(a & s ) diagrams 
corresponding to |(cc)8( 3 -Pj)g) vanishes for m q — > 0.) Since v 2 ~ 0.3 for charmonium, 
these contributions from the higher Fock states can easily be as large as (or even larger 
than) the electromagnetic decay of the valence Fock state. From these considerations it 
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should be clear that the standard extraction of the pion form factor from J/if> — > tttt is 
dubious. 

7 Summary 

In this paper we have presented a detailed analysis of Xcj decays into light pseudoscalar 
mesons within the framework of the modified HSA. For sufficiently heavy-quark masses, 
quarkonia are almost non-relativistic, and corrections to the quark-potential model de- 
scription can be organized into an expansion in v, the typical velocity of the charm quark 
in the meson. Recently it has been shown || that, for inclusive decays, the contribution 
from the higher Fock state |cc 8 ( 3 S'i)g) is not suppressed by v with respect to the contri- 
bution from the valence Fock state |cc"i ( 3 Pj)) (i.e. the quark-potential- model Fock state). 
Here we found a similar situation for exclusive Xcj decays: the colour-octet contribution 
(i.e. the one arising from the higher Fock state |ccg) with cc in a colour-octet state) is not 
suppressed by powers of either v or l/m c . Hence the usual suppression of higher Fock 
states for exclusive reactions does not hold in this case. 

Owing to the v -expansion, the wave function of the |ccg) Fock state of the Xcj is, to 
leading order in v, completely fixed apart from a single long-distance parameter, namely 
the colour-octet decay constant f( 8 \ However, in our approach J® remains the only 
parameter of the colour-octet contribution to the decay rate Xcj — ► MM (M = it, K, 77). 
We have given arguments against the use of higher Fock states for the pion. Therefore 
we conserve colour (an issue irrelevant for inclusive decays) by coupling the constituent 
gluon of the |ccg) Fock state of the Xcj to the hard process. This becomes a sensible 
proceeding in the mHSA. The two reasons for this are as follows. 

First, we keep the transverse momentum of the pion's valence-quark momenta. In this 
way, no singular integrals occur. The regularization depends on the pion wave function, 
but the latter is quite well constrained from an analysis of the photon-pion transition form 
factor in the mHSA. Secondly, incorporation of a Sudakov factor ensures that the coupling 
a s of the constituent gluon to the hard process never becomes large. The coupling is 
evaluated dynamically and the Sudakov factor suppresses soft phase-space regions where 
the scale of a s would become so small that a perturbative treatment were no longer 
justified. 

In our estimates of the decay rates Xcj - > KK, and 7777, the colour-octet con- 
tribution turned out to be of great importance in order to establish contact with the 
experimental data. The value of the single long-distance parameter came out very 
reasonable, consistent with the expectation from v scaling. Since the |QQg) Fock state of 
Xqj is neither v- nor power-suppressed, we expect also a large fraction of the Xbj — > 
widths to originate from the |bbg) state. This is indeed the case in our approach. Pre- 
dictions for bottomonium decays are quoted in Table 4. At present there is no data to 
compare with. 

In the case of exclusive J / ip decays, contributions from higher Fock states (and, hence, 
from colour-octet contributions) first start at 0(v 4 ) and are often also power-suppressed. 
They can therefore be neglected in most reactions, e.g. in baryon-antibaryon decay chan- 
nels, which are dominated by the contributions from cc annihilations through three gluons. 
Indeed, a recent calculation |31[] along the lines proposed here provides good results for 



many BB channels. We have pointed out that the situation is different for the J / ip decay 
into two pions. This decay is customarily assumed to be dominated by the electromag- 
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/( 8 )[GeV 2 ] 


T[Xbo -> 7r + 7r ] 


T[Xb2 ~> 7T + 7T ] 


T[Xbo - vr°7r°] 


T[Xb2 - vr°7r°] 


1.46 x 10- 3 


20.6 


1.69 


10.5 


0.88 


5 x 1(T 3 


152 


13.8 


78.2 


7.27 



Table 4: Decay widths in eV for P-wave bottomonia for two values of (m& = 4.5 GeV; 
P' p (0) = 0.7 GeV 5 / 2 ). 

netic decay of cc annihilation into a photon, since the three-gluon contributions cancel to 
zero. We have described the dominant hadronic decay channels arising from higher Fock 
states. These are likely to be responsible — at least partially — for the large difference 
between the value of the pion form factor in the space-like region and its value deduced 
from J/i[) decays into pions. Our considerations may also have consequences for the de- 
cays ip(nS) — > pn and for the process 77 — > im where, within the HSA, a substantial part 
of the cross-section is related to the pion form factor in the time-like region. 

A Appendix: Calculation of the colour-octet contri- 
bution within the mHSA 

The colour-octet decay amplitude can be written in a form similar to (3.4) : 

M^(xcj - Tr+Tr-) = jT 1 dxdy J K(y, ba) b x ) 

10 

x T l {J \x,y,b 1 ,b 2 ) exp [S(x, y, b x , b 2 , t a , U2)] • (A.l) 
i=i 

The sum runs over ten groups of graphs (see Figs. |4] and [|) in the case of charged pions. 
For the 7r°7r° final state the graphs of group 11 contribute as well. Each group contains 
a certain number of Feynman graphs, which differ only by permutations of either the 
two-pion states, the light quark and antiquark lines, or the attachment of the gluons from 
the charmonium side. Note that for those Feynman graphs where the cc pair annihilates 
into one gluon, the momentum fractions of c and c do not appear explicitly and hence 
only the momentum fraction z carried by the constituent gluon remains. For those graphs 
where cc annihilate into two gluons (groups 8-10) we assume that c and c carry each a 
momentum fraction (1 — z)/2 of the charmonium momentum. 

In Eq. (A.l) the T/ J ^ denote the Fourier-transformed hard scattering amplitude of 
group i for the colour-octet decay of the Xcj m transverse coordinate space. The hard 
scattering amplitudes T/ J ' ) are given in Figs. § and [|. 

In the expressions for the given in Figs. § and |5|, we absorb the combinatorial and 
colour factors as well as the colour-octet Xcj decay constant and a s into the factors 

64tt 5 / 2 / (8) I 

K = ~ 77? S— «sfe)«s(^2) V a s(*i3) , 

9V6m^ 

647T 5/2j(8) 

«2 = ^ — — a 8 (ta)a a (t i2 ) V^sfe) • (A.2) 



19 




Group 1 



2 - (l-z) 2 M 2 



g 2 = (z- Xl )(z- yi )M 2 



Ki 



9? 



z(z-xi)M 2 - 
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g 2 = {z-x 1 )(z-y 1 )M 2 -K 2 L 

q\ = z(z-x l )M 2 - k^j 



Group 2 
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Figure 4: Overview of groups 1-5 of Feynman graphs contributing to the colour-octet 
decay amplitude (M = 2m c , x\ = x, x 2 = 1 — x, is defined in (3.7)). 
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Figure 5: Overview of groups 6-11 of Feynman graphs contributing to the colour-octet 
decay amplitude (M = 2m c , x\ = x, x 2 = 1 — x, Kj_ is defined in (3.7)). 



21 



The renormalization scales appearing in the various amplitudes, depend on the kine- 
matics of the specific group. The virtualities of the two hard internal gluons connecting 
the external lines determine the renormalization scales tn and fa, respectively. In anal- 
ogy to ( p.ll|) these scales are chosen as the maximum of 1 /£>| and the corresponding 
gluon virtuality (at k_u = kj_ 2 = 0). In addition there is the constituent gluon from the 
Xcj colour-octet Fock state coupling to the hard part. The relevant scale t i3 is chosen as 
the maximum of 1/6?, l/frf. and the virtuality of the adjacent internal quark or gluon line 
(at k^i = k ± 2 = 0). 

We take into account the k± dependence only in the denominators of the T^ J \ i.e. 
we neglect corrections of 0(k"j_/m 2 ). In the collinear approximation we employed in M, 
propagator singularities are touched within the range of the momentum fraction integrals 
which cannot be regularized in the usual way with the te prescription. The difficulty arises 
from the gluon propagator g 5 (see Figs. f| and||). In || we regularized the corresponding 
integrals by replacing Ki in g\ through the r.m.s. transverse momentum (k^) 1 / 2 of the 
quarks inside the pions. Here, in the modified mHSA all integrals appearing in (A.l) 
are either finite or can be regularized with the ie prescription. In some of the T/ J ^ there 
are x-dependent terms in the numerators, which will cancel some of the propagator de- 
nominators if the transverse momentum is neglected. In these cases we actually let these 
terms cancel, thus avoiding some singularities. This procedure is equivalent to neglecting 
0{k\/m 2 ) terms. It is applied, for instance, in group 2 through which only the depen- 
dence on the transverse momentum difference K i (see (|T7P) is left in T 2 (J) . As already 
encountered in the case of the colour singlet case the corresponding Fourier-transformed 
hard amplitude T 2 ^ then includes a 5-function 5^{b\ — \>2). Similar situations occur in 
group 4 and in parts of groups 1, 3, 5 and 8. 

The Fourier transform of propagators regularizable with the te prescription read 

d 2 k ± expHk^-b] _ f-fH^yife) fors>0 
(2tt) 2 s-k 2 ± + te I i K (V=i6) fors<0 ' 



2tt 

Note that the inclusion of kj_ weakens the 1/s singularity to a logarithmic one. In the gen- 
eral case of two independent b- variables the integration over the relative angle 9 between 
hi and b 2 can be analytically performed by means Graf's theorem. For instance, 

/ dflH^tv/Jlbi-bal) = V / (Wcos(Z0)H, (1) (Vimax(6i,62))J J (v^niin(6 1 ,6 2 )) 

= 27rHf ) (Vs max(6 1 , b 2 )) J (v^ min(&i, b 2 )) . (A.4) 

In most of the cases one thus finds analytical expressions for the , except for group 1, 
where t[ j ^ depends on kj_i, k^ 2 and Kj_. Here, T[ J ^ includes the integral 



1 b db U^(y/(z-x)(z-y) + zeMb ) 



x W l, (Jz(z-x) +ieAfmax(6i,6o))Hj >(J(l-z)(z-y) + te M max(b 2l b )) 



x 3 (y/z(z-x) +zeMmin(6i,6o)) J (\/ (l-z)(z-y) + te M mm(b 2 ,b )) . (A.5) 

To proceed we split the range of integration at some large value bo of bo. Actually bo is 
chosen in such a way that all arguments of the Bessel functions Hq 1 '* and Jq are greater 
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than 2 for b Q > b , through which all Bessel functions can be replaced by their asymptotic 
expressions and (A. 5) is solved. For the remainder of the integration region (bo < bo) the 
integral is evaluated numerically. 
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